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m : ABSTRACT 
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i The structure of obscuring matter in the environment of active galactic nu- 



o 



clei with associated nuclear starbursts is investigated using 3-D hydrodynamical 
simulations. Simple analytical estimates suggest that the obscuring matter with 
^ I energy feedback from supernovae has a torus-like structure with a radius of sev- 

O I eral tens of parsecs and a scale height of ~ 10 pc. These estimates are confirmed 

by the fully non-linear numerical simulations, in which the multi-phase inhomoge- 
neous interstellar matter and its interaction with the supernovae are consistently 



O 

i-{ ■ followed. The globally stable, torus-like structure is highly inhomogeneous and 



6 ■ 

^ I turbulent. To achieve the high column densities (> 10^^ cm ^) as suggested by 

observations of some Seyfert 2 galaxies with nuclear starbursts, the viewing an- 
^ ■ gle should be larger than about 70° from the pole-on for a IO^Mq massive black 

' hole. Due to the inhomogeneous internal structure of the torus, the observed 

column density is sensitive to the line-of- sight, and it fluctuates by a factor of 
order ~ 100. The covering fraction for N > 10^^ cm~^ is about 0.4. The average 
accretion rate toward i? < 1 pc is ~ O.AMq yr^^, which is boosted to twice that 
in the model without the energy feedback. 

Subject headings: ISM: structure, kinematics and dynamics — galaxies: nuclei, 
starburst — method: numerical 
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1. INTRODUCTION 

Optically thick obscuring molecular tori have been postulated to explain various proper- 
ties of active galactic nuclei (AGN), especially the two major categories of the AGN, namely 
type 1 and type 2. However, the true structure and the formation and maintenance mecha- 
nism of the torus have not yet been understood either theoretically or observationally. 

Recently a number of observations have suggested a new aspect of Seyfert 2 galaxies 
(Sy2). It is pointed that there are possibly two types of Sy2s: one is a classic Sy2, an obscured 
Seyfert 1 nucleus. The other type II is a Sy2 with a nuclear starburst (Cid Fernandes & 
Terlevich 1995; Granato, Danese, & Pranceschini 1997; Gonzalez Delgado, Heckman, & 
Leitherer 2001; Levenson, Weaver, & Heckman 2001). In the latter case it is believed that a 
nuclear starburst is associated with the obscuring material whose scale is i? < 100 pc. These 
observations imply that the classic picture of the unified model for AGNs, in which the 
diversity of the Seyfert galaxies is explained only by the geometrical effect of a geometrically 
and optically thick torus, may not be applicable to some fraction of the Sy2s (see also 
Maiolino et al. (1995)). More generally we propose here that there is an additional parameter 
in the unified models, namely, the strength of the nuclear starburst. It is this nuclear 
starburst and the mass of the black hole (BH) that determines the geometry of the obscuring 
torus. 

While discussing the X-ray background, Fabian et al. (1998) proposed that low- luminosity 
AGN could be obscured by nuclear starbursts within the inner ~ 100 pc. They suggested 
that supernovae (SNe) from a nuclear starburst can provide the energy to boost the scale 
height of circumnuclear clouds and so obscure the nucleus. For simplicity, they assumed 
that the material sits in an isothermal sphere, and the gravitational effects of the central 
BH are ignored. This idea seems to be consistent with the recent observations of Sy2s with 
nuclear starbursts. Here we investigate the model in more detail. The problem is actually 
very complicated. We need to solve the 3-D inhomogeneous structure and dynamics of the 
ISM in the combined gravitational potential of the central massive BH, the central stellar 
system and the central massive gas distribution. We must handle the radiative cooling not 
only for the hot gas, but also for the cold gas, and include the feedback interaction between 
SNe and the inhomogeneous ISM. Numerical simulations are powerful tools for this kind of 
complicated problem. 

Recently we presented high-resolution, 3-D hydrodynamical modeling of the ISM in the 
central hundred pc region in galaxies, taking into account self-gravity of the gas, radiative 
coofing, and energy feedback from SNe (Wada & Norman 2001; Wada 2001) (hereafter WNOl 
and WOl). In this Letter, we apply this numerical method to the gas dynamics around the 
central massive BH, and study the starburst- AGN connection in Seyfert galaxies. 
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2. GEOMETRY OF OBSCURING MATERIAL AROUND THE CENTRAL 

ENGINE 

In this section, we discuss the possible structure of the obscuring material around the 
central massive BH under the influence of the SN explosions detonating in the central mas- 
sive gas cloud. The obscuring material will be inhomogeneous, and multi-phase and the 
geometrical thickness is assumed to be supported by its internal turbulent motion caused by 
energy input from SNe. 

We assume the turbulent energy dissipation in a disk around the central massive BH 
is in equilibrium with the energy input from SNe. In a unit time and volume, the energy 
balance is 

^ = T^s^E, = m-PgEs. (1) 

where Vt is turbulent velocity of the gas, is the dissipational time scale of the turbulence. 
Eg is the total energy injected by a SN, 77 is a heating efficiency per unit mass that represents 
how much energy from SNe is converted to kinetic energy of the ISM, and Sy^ is star formation 
rate per unit volume and time. We can assume the star formation rate is proportional to 

the dynamical time scale in the disk namely ~ aVtcTjgh~^ = aVc/rpg, where Vc is the 
circular velocity, r is a radius, h is the scale height of the disk, is surface density, and 
a is a constant (a ~ 0.02 is suggested for a global star formation in spiral galaxies rate 
by Kennicutt 1998). The dissipational time scale can be simply assumed as = ^h/vt, 
where ^ is a constant of order unity. Therefore eq.(l) is 

I = Va'-Eg. (2) 

Assuming hydrostatic equilibrium for the vertical direction of the disk, i.e. pgV^ ~ Pggh, we 
have ~ G'Mbh'" for the case that the central massive BH dominates the gravitational 
potential (regime I), or f| ~ ttGS^/i, if the stellar disk or bulge potential dominates the 
potential. Here is surface mass density of the stellar system (regime II) . For the regime 
I, h/r ~ Vt/vc and for the regime II, h/r ~ Vt/'^c- Therefore from eq.(2), we have h ~ 
{^rjaEsY^'^rv'^, for the regime I, or for the regime II, h ~ {^'qaEsYrv''^. Using Vc,i — 
(GMbhA)^^^ for the regime I (r < ro), or Vc,ii = (7rG'Ey,)^/^r^/^ for the regime II (r > ro), 
the disk scale height in each regime is 

hj{r) = (e77a£;.)V^G'-V2M-^/V/^ r < ro (3) 

and 

hii{r) = (e?7«£;,)2(7rG'E,)-V-\ r > ro (4) 
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where ro is determined from Vcjiro) — Vc,ii{ro), and it is ro = 

where the BH mass. Ms = lO^M© and E*_4 = E*/lO^M0pc~^. Therefore the maximum scale 
height at each regime is 

^ 22pc {i,a,^^,E,,fl^Ml'\T,T.,,,)-'l\ (5) 

^ lpc(6a2r;-3^5i)X"'^'s;f , (6) 

where 6 = C/1-0, 0:2 = q;/0.02, 77_3 = 10-^M-\ the total gas mass Mg^s = Mg/IO^M^. The 
scale height can be also expressed as /io,/(^o) ^ 35 pc S F R.y^rl^'^ {^lE^i/ Mg^^y^'^ , with the 
star formation rate SFRi = lAf© yr~^ and rg = r/60 pc. Here we assume the Salpeter's 
IMF. Note that /iq,// > Kj, if ^2 > 0.08 for ^1, r/_3, Mg, S^,4, and E^i. Eqns (3) and (4) 
show that the disk has a torus-like structure (it looks more like a concave lens rather than 
a bagel), and its opening angle $ is tan$ ~ ro/ho,i oc (,^?7a£'5)~^/^Mg^ , for ho,ii < /io,/- 
Therefore $ ~ 70° for ^1, 77-3, 02, Afg, and S^^4. For the case that nuclear star burst is very 
active, e.g. SFRi — 10 and 100, we have smaller opening angles, i.e. $ ~ 50° and 20°, 
respectively. 

If the star formation originates via gravitational instability in the gaseous disk, stars 
are not formed near the massive black hole. Assuming the Toomre's stability criterion 
'^crit — ^cCs/t^G — CsG~^/'^M^^r~^/'^ /n, where Cg is the sound velocity, the critical radius Tc 
for which the disk with a uniform surface density is unstable is Tc ~ [c^rgMBH/ {MgG)\ , 
then Tc ~ 3pc c^^r^^ mI^^ Mg^^^ , where ci = Cs/{1 km s~^). The disk scale height for 
the stable disk would be very thin inside r — rc, otherwise its shape would be represented 
by eqns (3) and (4). Therefore the gas around the central massive black hole with the 
nuclear starburst naturally forms a geometrically thick obscuring structure around the AGN, 
depending on the mass of the black hole, the gas mass, and the strength of the burst. 

For the case of hoj ^ hoji, the average column density toward the nucleus as a function 
of the viewing angle (0 = is pole-on), N{(f)) ~ {pg)L{(j)), where L{(j)) is the path length in 
the torus, is approximately N((f)) ~ 5x lO^^cm"^ Mg^s^'^^h22^{(j))., where /i22 = h/22 pc. r(0) 
is a function that determines the path length in the torus, r(0) ~ (1 — 2tan~"'^ 0)"*^^^/ sin0, 
for To/hoj < tan0 < ro//io,//, or r(0) ~ cos""*^ (p, for tan0 > ro/hQjj. F ~ 0.8 for (p ~ 80°. 
If the star formation rate is ten times larger and hoj < hojj, N{(f)) is smaller by a factor of 
two. 
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3. NUMERICAL MODELING OF THE OBSCURING GAS IN THE 

GALACTIC CENTER 

The argument in §2 is basically an order- of- magnitude estimate, and, consequently, 
detailed structure and evolution of the ISM in the galactic center should be verified by 
a fully non-linear, time-dependent numerical simulations. The star-forming massive gas 
should be an inhomogeneous with a velocity field that is turbulent (see WNOl and WOl). 
SN explosions in such inhomogeneous system cannot be simply parameterized. The evolution 
of the SNe remnants will not be spherically symmetric, and should not be expressed by a 
simple analytic solution. As a result, it is not straightforward to estimate how much energy 
of a SN explosion is converted to the kinetic energy of the inhomogeneous, non-stational 
ISM. The energy dissipation rate as a function of time and position in the non-stationary 
turbulent gas disk needs to be explicitly calculated. In this section, we use a state-of-the-art 
numerical technique to clarify the evolution and structure of the gas in the star-forming 
region around the central massive BH. 

3.1. Numerical Method and Model 

The numerical methods are the same as those described in WNOl and WOl. Here we 
briefly summarize them. We solve the same equations as eqns.(l)-(4) in WOl numerically 
in 3-D to simulate the evolution of a rotating ISM in a fixed gravitational potential. Here the 
external potential force in the momentum conservation equation (eqn. (2) in WOl) is V$ext + 
V$BH, where the time-independent external potential is $ext = — (27/4)^/^'t;^/(r^ + a^)^/^ 
with the core radius a = 10 pc and the maximum rotational velocity Vc — 100 km s~^. The 
central BH potential is $bh = —GMBu/{r'^ + b'^Y^'^ with 6=1 pc. We also assume a cooling 
function A{Tg) (20 < Ty < lO^K) with Solar metallicity and a heating due to photoelectric 
heating, Fyv and due to energy feedback from SNe, Fsn. We assume a uniform UV radiation 
field, which is ten times larger than the local UV field. The hydrodynamic part of the basic 
equations is solved by AUSM (Advection Upstream Splitting Method) (Liou & Stcffcn 1993). 
We use 256^ x 128 Cartesian grid points covering a 64^ x 32 pc^ region around the galactic 
center (the spatial resolution is 0.25 pc). The Poission equation is solved using the FFT 
and the convolution method. The initial condition is an axisymmetric and rotationally 
supported thin disk with a uniform density profile (thickness is 2.5 pc) and a total gas mass 
of Mg — 5 X IO^Mq Random density and temperature fiuctuations, which are less than 1 % 
of the unperturbed values, are added to the initial disk. 

SN explosions are assumed to occur at random positions in the region of |a;|, |y| < 51 pc 
and |-z| < 4 pc. The average SN rate is ~ 1 yr~^, which is corresponds to a ~ 0.3 at r = 20 
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pc. The energy of 10^^ ergs is instantaneously injected into a single cell as thermal energy. 
The 3D evolution of blast waves caused by SNe in an inhomogeneous and non-stationary 
medium with global rotation is followed explicitly, taking into account the radiative cooling. 
Therefore the evolution of the SNRs, e.g. the duration and structure of the SNe, depends 
on the gas density distribution around the SN. 

3.2. Numerical Results 

Figure 1 (a), (b) and (c) show volume-rendering representations of the density distri- 
bution of the gas around the massive BH in a quasi-stable state (t = 1.6 Myr) from three 
different viewing angles. In Fig. 1(a), the nuclear region can be seen, whereas the nucleus is 
obscured in Fig. 1(b) and 1(c). The right half of the Fig 1 (b) is a cross-section of the thick 
disk, and the inner structure of the thick disk is clear; the scale height of the disk is smaller 
in the central region than in the outer region as expected in §2, and the density structure 
is very inhomogeneous. The three panels in Fig. 1 together show that the gases around the 
central massive BH form a torus-like structure. The inhomogeneous, filamentary structure 
is quasi-stable. That is, the local inhomogeneous structure is time-dependent, but the global 
torus-like morphology does not significantly change during several rotational period. 

In Fig. 2, we plot line-of-sight column density of the gas as a function of the viewing 
angle (90° is edge-on). The column density is higher in edge-on view. To achieve N > 10^^~^^ 
cm^^, which is suggested in observations of Sy2s with nuclear starbursts (Levenson, Weaver, 
& Heckman 2001), the viewing angle should be larger than ~ 70° in this model. It should 
be noted, however, that the column density has large (about two orders of magnitude) 
fluctuations for the same viewing angle. This is caused by the inhomogeneous internal 
structure of the torus, which shows a Log- Normal distribution function (see WOl). A factor 
of 100 change in the average column density corresponds to about ±20° in the viewing angle. 
We find that a fraction of the sohd angle, for which the nucleus is obscured, is 0.21, 0.37, 
0.46, and 0.55 for N > 10^^, 10^^, 10^^, and 10^^ cm~^, respectively. These numbers are a 
factor 2-3 smaller to explain the ratio of Sy2s to Syls, which is about 4 (Maiohno & Rieke 
1995), based on the 'strict' unified model. 

From eqns. (5), we can expect the scale height of the disk of for a = 0.3 is /i ~ 16 pc 
at r = 20 pc and /i ~ 6 pc at at r = 10 pc. This consistent with the numerical result. The 
column density predicted in §2 for (f) ~ 80° fit the simulation. However, Fig. 2 shows that the 
column density depends more strongly on the viewing angle than the analytical prediction. 
This is because the average density distribution for a ^-direction in the torus is not constant, 
contrary to the assumption in §2. 
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Fig. 3 shows the time evolution of the gas mass inside r = 1 pc for the models with and 
without energy feedback. The average accretion rate is Mg ~ OAMq yr~^, which is roughly 
twice that in the model without feedback. The turbulent viscous time scale in this system 
is ttub ~ f^^/'^tuh with viscous coefficient z/tub- Using z/tub ~ "^t x h, where If is the largest 
eddy size which is about the scale height of the disk (~ 10 pc), the turbulent viscous time 
scale is ttub ~ 10^(t't/10kms~^)~'^ yr. This is the same order of the numerical result, i.e. 
Mg/Mg - 5 X lO^M0/O.4M0yr-^ 

4. DISCUSSION 

The estimate in §2 suggested that the opening angle positively correlates with Mbh- 
As Mbh increases, the critical radius for the gravitational instability is also larger. This 
implies that low luminosity (i.e. small Mbh) AGNs are obscured for larger fraction of the 
viewing angle than the high luminosity AGN, if, quite reasonably, the mass accretion rate 
has a positive dependence on the mass of the central engine. In other words, luminous type 
2 objects, such as type 2 quasars are less frequently observed in the optical than type 1 
quasars. They could be observed mainly X-rays, and this is consistent with observations 
(Norman et al. 2001; Lehmann et al. 2001; Akiyama et al. 2000). 

Krohk & Begelman (1988) (KB) claimed that stirring by stellar processes is never strong 
enough to compete with energy dissipation in the clumpy torus. One should note, however, 
that in their order of magnitude estimate, a much smaller, thicker torus (~ 1 pc) is assumed 
than in our extended disk model. Therefore KB need large velocity dispersions (> 200 km 
s~^) to keep the disk thick, because thick disks with h/r ~ 1 require Vt/v^ ~ 1. The velocity 
dispersion for /i/r ~ 1 can be less than 50 km s^^ for r > 30 pc. The smaller velocity 
dispersion is favorable, because the energy dissipation rate is smaller (see eq.(2)). KB also 
pointed out that a luminosity-to- gas-mass ratio, L/Mg, for the stellar stirring model would 
be a factor 100 larger than the observed values. This is also not the case for a 100 pc 
scale obscuring material. For the analytical model presented here, L/Mg ~ r~^/^. In fact, 
a much larger molecular gas mass is expected on such a scale (e.g. Kohno et al. (1999)). 
Recent IR and X-ray observations of Scyfert nuclei (e.g. Granato, Danese, & Franceschini 
(1997)), on the other hand, suggest a 100 pc scale extended obscuring material. This picture 
is consistent with our model of the extended starburst supported obscuring region in which 
case we still obtain a geometrically thick disk with velocity dispersions < 50 km s^-*^. Note 
also Pier & Krolik (1992) and Ohsuga & Umemura (2001) for the effects of radiation pressure 
on supporting the obscuring torus. 

The present model suggests that the ratio of Sy2s to Syls cannot be explained only by 
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the orientation of the viewer relative to the torus. In Syls, the star formation rate in the 
nuclear region would be smaller than in Sy2s. The strict unified model, in this sense, should 
be modified by including the strength of the nuclear starburst. 

We thank our colleagues Roberto Cid- Fernandez, Tim Heckman, Julian Krolik and 
Nancy Levenson for useful and stimulating discussions. Numerical computations were carried 
out on Fujitsu VPP5000 at NAOJ. KW is supported by Grant-in- Aids for Scientific Research 
(no. 12740128) of JSPS. 
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Fig. 1. — 3-D density distribution represented by a volume-rendering technique. The colors 
represent the relative opacity based on the line-of-sight column density, but it does not 
represent the absolute opacity or any physical value. The regions colored red or yellow are 
more optically thick than the blue regions. The right half of (b) is a surface section. 

Fig. 2. — Column density vs. viewing angle. The column density toward the nucleus is 
calculated for the density distribution of a model shown in Fig. 1 every ^ ~ 4° and 4° for 
azumithal direction. 

Fig. 3. — Time evolution of the gas mass inside i? < 1 pc for two models (with and without 
energy feedback). The solid hue represents the mass accretion rate O.3M0 yr"^. 
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